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It is shown that a nonequilibrium environment can be instrumental in suppressing decoherence 
between distinct decoherence free subspaces in quantum registers. The effect is found in the frame- 
work of exact coherent-product solutions for model registers decohering in a bath of degenerate 
harmònic modes, through couplings linear in bath coordinates. These solutions represent a natural 
nonequilibrium extension of the Standard solution for a decoupled initial register state and a ther- 
mal environment. Under appropriate conditions, the corresponding reduced register distribution can 
propagate in an unperturbed manner, even in the presence of entanglement between states belonging 
£\J \ to distinct decoherence free subspaces, and despite persistent bath entanglement. As a byproduct, 

we also obtain a refined picture of coherence dynamics under bang-bang decoherence control. In 
' particular, it is shown that each radio- frequency pulse in a typical bang-bang cycle induces a revival 

of coherence, and that these revivals are exploited in a natural way by the time-symmetrized version 
of the bang-bang protocol. 
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í> . I. INTRODUCTION 

At a very fundamcntal level the control of decoherence in a quantum register amounts to the engineering of its 
£SJ | entanglement with the surrounding environment. The earliest approach to the problem was guided by the philosophy 
of classical error correction codes, and has culminated into a comprehensive theory of quantum error correction codes 
[QECCs] In essence QECCs are designed to actively purge unwanted environment entanglement by a recurring 

error correction cycle, to the cffcct of recovering periodically an 'error-free' reduced register state. Other active 
strategies, known as quantum bang-bang control or dynamical decoupling j^-^J], have been developed by analogy with 
Oh, NMR decoupling and refocusing techniques. In this case entanglement with the environment need not necessarily 
■ cancel at the end of the correction cycle, but is coerced, by suitably engineered control "pulses", to "average out" 
in such a way as to leave the reduced register state virtually unchanged p5[ . Alternatively, it was shown that 
specific forms of decoherence can be passively suppressed, provided the intrinsic symmetries of the noise-generating 
interactions allow the existence in the register Hilbert space of decoherence- free subspaces [DFSs] that are inaccesible 
to noise processes jï^-|20|. Since any entanglement with the environment leaves a DFS code unaffected, robust storage 
of the reduced register state can be achieved without further active correction cycles. The concept of DFS has emerged 
recently as only a particular instance of symmetry generated noiseless subsystem [NS] structures, capable to support 
robust passive storage through redundant encodin g |2l| 1 . Both types of structures have been confirmed experimentally, 
through the implementation of optical DF states [|22||, of one-qubit DFS memòries in two trapped ions p3], or in two 
nuclear spins p3[ |, and of a one-qubit NS memory encoded in three nuclear spins pEfl . Notably, QECCs and dynamical 
decoupling strategies have been also shown to enforce symmetrized effective evolutions capable to support NSs, so 
that noise-protected encoding takes place in fact in dynamically generated NSs [^6| . From a broader perspective, in 
all cases decoherence control aims to maintain, or recover periodically, a given reduced register state in the presence 
of time-dependent entanglement with environmental modes. The ability to implement such effective noise protection 
strategies is paramount for the practical development of both robust quantum memòries and fault-tolerant quantum 
computation. 

For a typical model register of N two-level elements, immersed in a quantized environment and described by 

JV-l N-l 

H =J2 e ^ n) + E nLj i b l b i + EE ^ hc qn b\ + x* qn h) > (i) 
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the passive or active character of the control strategy seems to be conditioned by the number of DFSs accessed by the 
encoding method. Indeed, since the associated Hilbert space is always decomposable into a direct sum of elementary 
DFSs [at the very least, direct products of the 1-dimensional eigenspaces for each o* ], the corresponding reduced 
register state is necessarily distributed over a [finite] number of disjoint DFSs. If the encoding state is confined to a 
single DFS, decoherence is passively suppressed, even under entanglement with the environmcnt. But if the encoding 
involvcs múltiple entangled DFSs, active control seems mandatory. 

The main intention of this paper is to point out that decoherence-free states, and therefore passive control, may be 
possible as well for nonseparable distributions entangled over múltiple DFS, and also entangled with the environmcnt. 
Such distributions realize a noise-protected concatenation of distinct DFS, which suggests, in principle, that the 
capacity for coherent storing of a DFS-supporting register may extend beyond the dimensionality of the largest DFS. 
Sec. II presents a general argument in favor of this possibility. For a more precise discussion, we exploit an extended 
class of exact, closed-form density matrix solutions for systems with decoherence-free subspaces under interactions 
linear in environmental degrees of freedom, which have been retrieved as density matrix generalizations of the pure 
state coherent-product [Davydov] ansatz in ref. |^7j. After providing some straightforward Davydov [pure] state 
examples of decoherence-free propagation on múltiple DFS, in Sec. III we introduce the thermal coherent-product 
states in a self-contained manner, and with reference to the specific context of quantum registers. These states 
involve a nonequilibrium environment in a statistical superposition of Gaussian states and provide a natural extension 
of the Standard solution currently employed in discussions of the decoherence process |^,^], which refers to an 
uncorrelated thermal environmcnt. We show in Sec. IV that the corresponding register reduced density matrix can 
propagate, under proper conditions, in an unperturbed, decoherence free manner, even when nontrivially distributed 
over múltiple, nonisomorphic DFS, and strongly entangled with the environment. In Sec.V we discuss the conditions 
for multi-DFS decoherence free evolution for three particular cases of the spin-boson model ([!]): a single-qubit register, 
a multi-qubit register with weak collective decoherence and a multi-qubit register with individual decoherence. In 
Sec. VI it is pointed out that sudden spin-flips of the type considered in bang-bang protocols [0 imply in effect a 
manipulation of the register-environment entanglement that may be of interest for the preparation of multi-DFS 
decoherence free states. Alternatively, it is also found that the density matrix coherent-product solutions for model 
(Q) prové instrumental for an analysis of the detailed coherence dynamics under bang-bang control S . It is shown 
that each qubit spin-flip in a bang-bang cycle produces a revival of coherence, and that these revivals provide an 
alternative justification for the well-known superior performance of the time-symmetric version of the protocol [ p"2[ . 
A summary and conclusion are provided in Sec. VII. 



II. THE GENERAL MODEL 



Consider a quantum register R, described by a Hamiltonian Hr, and immersed in an environment B with Hamilto- 

nian Hr- Let the register-environment interaction be given as H$b — 8> Bk, with Rk [Bk] self-adjoint register 

fe 

[environment] coupling operators, such that the total Hamiltonian reads 

H = H R + H B + R * ® B * ■ ( 2 ) 

fe 

Let us assume also that the register dynamics supports múltiple decoherence-free subspaces Sj labeled by an appro- 
priate index set { J}. Here we refer to the closed-system, hamiltonian version of the theory of DFSs fïl||l9|| . That is, 
let A denote the f-closed register interaction àlgebra generated by the hamiltonian Hr, the register components of the 
interaction Rk, and the identity Ir on the register Hilbert space. By a bàsic theorem of f-closed operator algebras, 
A decomposes into a direct sum of dj x dj-dimensional irreps, each of multiplicity rij, i.e. 

- 4 = / «./ ®M(dj x dj,C) , (3a) 
j 

where A4(dj x dj, C) denotes the àlgebra of dj x dj complex matrices. As a result, the register Hilbert space Hr 
also decomposes into a direct sum of subspaces invariant under A, as 

Hr = € nj ® € dj . (3b) 

,7 

The register DFSs Sj correspond then to the 1-dimensional irreps JM(1 x 1, C) of A, and their dimension is identical 
to the multiplicity of the associated irrep, dimSj = nj. Equivalently, a necessary and sufficient condition for the 
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occurrence of DFSs pffl is that the kernels of all commutators [Hr, Rk] have a nontrivial intersection [e.g., [Hr, Rk] = 
0, as in model (0)]. Further, 

i) any register state |?/>j ) belonging to the J-th DFS, IV'j ) € Sj, is a simultaneous eigenstate of all register couplings 
Rk for real eigenvalues /ifc.j common to the entire DFS, 

R k \^) = nkMj) ; 

ii) the DFS is left invariant by the register hamiltonian Hr, 

H a \1>f) € Sj . 

A ubiquitous example of register supporting DFSs is provided by a spin-boson model (üh with weak collective deco- 
herence, i.e., with identical coupling of each spin to a common bath, via [e.g.] the eri component only, such that 
Xqn = Xq [the morc general collective decoherence case involves all spin components, through interactions of the form 



?(") . 



N—í 



Xq\ 



In this case the register-bath interaction reduces to S Z Y^ {Xqb\ + Xq^q), where S z = a z denotes the 

q n=0 

z-component of the total spin, and the DFSs are the [highly] degenerate eigenspaces of S z . 

Returning now to our general model, throughout the following we shall restrict the register state space to 

7e = 05j , 

./ 

such that if B is the state space of the environment, the joint register-environment states live in the direct product 

Hilbert space Tí = 1Z <£> B. In this case, the subalgebra of A generated by the restrictions of Hr, Rk and Ir on 1Z 

is necessarily abelian. Hence, if Pj are the self-adjoint, orthogonal projectors on the DF subspaces Sj, such that 

PjPk — PrPj — òj t K, and Pj = ^TCj where I-r is now the identity on 1Z, properties (i-ii) above show that 

J 

i') every Rk has on 1Z the simple form 



R k 



j 



Vk,jPj , 



(4) 



ii') [Hr,Pj] = and, accordingly, the total Hamiltonian (||) describes a quasi-separable dynamics on 1Z. That is, 
it generates a superposition of independent evolutions on individual DFSs, where each DFS creates its own separable 
environment. Indeed, Hamiltonian (Q) can be written, by making use of the orthogonality and completeness of the 
projectors {Pj} on 1Z, 



H 



E 

j 



PjHrPj + Pj ® [h b + Ç Vk,jB k j 



and the evolution it generates obviously factores into separate evolutions on the individual DFS. 

Our starting observation is that, alternatively, the entire unperturbed register evolution can be formally factored 
out of the total evolution. It is sufficient to note that the register Hamiltonian Hr = ^ PjHjPj commutes with each 

J 

of the environment driven terms in the sum. Then the evolution operator generated by H can be conveniently cast as 













exp 




= exp 




■ exp 



H 



B,J 



exp 



-n HRt 



exp 

j 



-j.H BíJ t 



Pj , 



with 



Hb,j — Hb + 2^ Vk.jBk 



(5) 



(0) 



This factorization hints, in a rather evident manner, to the possibility of unperturbed and decoherence free register 
dynamics on múltiple DFS, since the register reduced density matrix at time t acquires the form 
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5*(t) 



exp 



-H R t 



■p (t) ■ exp 



■Hat 



with 



T{t)=T n 



exp 



~rH B ,jt 

n 



p(0) ■ exp 



n 



(7a) 



(7b) 



When the initial total density matrix p(Q) is such that at all times 

f(t)=p R (0) , 

the register evolution becomes indeed unperturbed and decoherence free, 



p R (t) = exp 



p R {0) ■ exp 



H HRt 



However, there is no indication yet that the conditions for this peculiar dynamics can be met outside of individual 
DFS. For a better picture, let the initial density matrix p(0) be given a representation in terms of a tensor product 
basis in Tí, of the form {\(f>%) <£> \(t>&)}, where {\4>a)} is an orthonormal basis in B and {\(f>%)} an orthonormal register 
basis yielding an irreducible representation [irrep] on 1Z, with v an index labeling basis states belonging to the same 
DFS. That is, let 



tf°) = E E i^>®i^)p^w(ox^i®(^i • 

Jv,(7 J'l/' \<j' 

Then the corresponding reduced density matrix for the register [Eq.(Qb)] can be brought to the more transparent form 



p H {t) = exp 







:Ï Hr \ 





Jv J'v' 







■ exp 


Ï Hr \ 



where the environment driven correlation amplitudes Xj v ,j' v '{t) are given by 

^Ju,.J'u'{t) = ^2 PJv,<r;J'v',cr'(Q)(4>f ' 



|exp 









h 


■ exp 


n 



\ti 



(8a) 



(8b) 



Expression (|8|a) now shows in fair detail that the reduced register evolution can become unperturbed, even in the 
presence of environment-mediated entanglement between distinct DFS [J ^ J'J, provided the correlation amplitudes 
^jv.j'v'it) remain constant in time. Note that any correlations \j v ,Jv> between register states in the same DFS 
are always constant in time, and the projected dynamics Pjpn(t)Pj on any DFS is efectively decoupled from the 
environment. This is a direct expression of the quasi-separable nature of Hamiltonian (^|). 

It is also worth pointing out that, by construction, such multi-DFS decoherence free states are distinct from dis- 
tributions on noiseless subsystems f2lf| . The latter are defined by nonabelian, multidimcnsional irreps M{dj x dj,C) 
of the register interaction àlgebra A [see Eq.(||a)] and, according to decomposition (||b), reside in the subspaces 
C™ J ® C dj , dj > 1, of the total register Hilbert space. Note that a DFS is in effect a trivial noiseless subsystem, 
supported by a direct sum of equivalent [identical] 1-dimensional irreps of the interaction àlgebra. In contrast, the 
multi-DFS decoherence free states referred to above are defined on a direct sum of prcsumably non-isomorphic DFS, 
i.e., on a direct sum of non-identical 1-dimensional irreps of the register interaction àlgebra. The register [subjspace 
1Z carrying these states is in fact orthogonal to any subspace corresponding to a nontrivial noiseless subsystem. 



III. THE COHERENT PRODUCT STATES 



A. The model environment 



From the expression of the density matrix factor (^>), or equivalently, from the correlation amplitudes (^b), it may 
be inferred that the exact nature of cross-DFS decoherence free distributions, if any exist, will depend considerably on 
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the particular environment interaction àlgebra generated by Hb, the couplings and the identity Ib on B. Foregoing 
a comprehensive analysis, we aim to prové that such distributions indeed exist in the familiar setting of a harmònic 
[bosonic] environment with linear couplings. Then Hamiltonian (|J) becomes, in the usual notation, 



= H R + J2 foqbpq +J2 Pj ® i^Jqbl + (J.*j q b q ) . 

i J,q 

where in the last line above the coupling constants have been redefined for convenience as 

l·Jq = /^kjXkq ■ 
k 



(9a) 



(9b) 



In this case, each DFS creates a decoupled environment of displaced harmònic modes, driven by Hamiltonians of the 
form 



H B ,J = hLJ q b l b q + (V>jq b \ + Vjq b q) = ^ ^^q 3 ^ ~ ' 



q q 
with statically displaced modes described by 



bj q =b q + 



Tllúa 



and energy shifts given by 



M 2 

hlÜn 



(10a) 



(10b) 



(10c) 



Accordingly, decomposition (ph for the total propagator becomes 



e ~i Ht = e~^ HRt Y e% n ° jt exp 



i 

h 



J2 nUJ 1 h lq h Jl 1 



(11) 



B. Simple examples of multi-DFS decoherence free states: pure coherent-product states 



Given the structure of the propagator ( |l·l| ) , it is not difhcult to identify pure states of the joint register-environmcnt 
system that generate nontrivial, multi-DFS decoherence free solutions. Consider, for instance, a total pure state |<I>) 
formed as a superposition of unnormalized [but orthogonal] register states |</>j) defined on [some] individual DFS, 
each correlated with the displaced vacuum |0j) of the corresponding environment Hamiltonian Hb,j [Eq.([ïo|a)1, i.e., 



(12) 



Its time evolution follows straightforwardly by applying expression ( |l·l"| ) for the total propagator, and reads 

I*) = Y e ~ ÍHRteÍn ° t ^ (g> l° J > ■ ( 13 ) 



Obviously, the associated register reduced state, given by the reduced density matrix 



p R (t) 



; *(^-^)*|^)( 0j ,| 0j )(^| 



(14) 
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evolves unitarily. Moreover, it evolves in an unperturbed, decoherence free manner, provided the energy shifts for the 
DFS involved are identical, that is 



2-^ hu>„ 



E 



\V£q\ 
TlLJ„ 



(15a) 



It also involves nonzero correlations between register states belonging to distinct DFS if the displaced environment 
vacua associated to those DFS have nonvanishing overlap, such that 



(0j'\0j) ^° 



(15b) 



One may recali that reciprocally displaced harmònic vacua generate unitarily equivalent Fock spaces iff their overlap 
amplitude is nonvanishing. Hence condition (|Ï5|b) requires in effect that the corresponding environment vacua must 
reside in unitarily equivalent spaces. That is, e.g., the Fock basis states |{ïïj'}) generated from \0j>) must reside in 
the Fock space spanned by the states generated from {|0j)}, and can be related to the latter by a unitary 

transformation that is well-defined on the subtended space, as = U(bj q , bj )\{nj}) [in particular, |0j<) = 

U(bj q , bj q )\0j)]. If the vacua are unitarily inequivalent, such that (0j<|0j) = 0, the reduced state (Q) does evolve 
unitarily, but as a trivial, separable superposition of distributions belonging to orthogonal DFS. 

A more sophisticated example is obtained if the DFS vacua \0j) are replaced in the total state (|ï^ ) by arbitrary 
coherent states of the corresponding displaced modes, defined by 

b Jq \{flj, 9 }) = Pj, 9 \{0j, 9 }) • 

In addition, let us allow for the possibility that every orthonormal register basis state \4>%) be associated with a 
different coherent state. The result is a Davydov coherent-product state 



Jv 



which propagates as 



|$) =X)^e-* ff «*e* n ° í |^,) ® |{/W*)}> 



(16) 



(17) 



Jv 



The states \{f3j v ,q(t)}) are coherent states evolved under the DFS-specific environment Hamiltonians Hb,j, and are 
dctcrmhicd, as usual, by 



0Jv,q(t) — Pjv,t 



The reduced register state reads now 



p H (t) = e~* HRt 



E ei(0 ° , " í2 °'' )t ^c} v ,|^)(/3 J v', ? (í)l/3^W>(^v'l 



Jv.J' v' 



■ pJ Hrí 



(18) 



and the corresponding conditions for nontrivial decoherence free propagation amount to 

i) the energy-shift condition (|Ï5|a); 

ii) the requirement that the interference amplitudes for coherent states belonging to different DFS be stationary in 
time, 



iii) the requirement that [some] coherent states \(3j Uiq (t)} associated to different DFS reside in unitarily equivalent 
Fock spaces, such that 

(/W >a (t)|/W*)> ^ • 

The second condition becomes more intuitive if recalled that the coherent states \(íj v ,q{t)) and \f3j' V >,q{t)) move on 
"orbits" centered on the reciprocally displaced vacua |0j) and |0j'). Assuming condition (ii) satisfied, condition (iii) 



G 



requires then that this movement be synchronized or phased so that the overlap (/3jv ,g(í)|A/i/,g(í)) remains constant 
in time. Note that coherent states associated to the same DFS are always synchronized and their overlap, if nonvan- 
ishing, is necessarily time- independent. 

Although an analysis of conditions (ii) and (iii) above is straightforward, we postpone following this direction. 
Instead, we proceed to construct a class of finite-temperature generalizations of these Davydov states. Subsequently 
it is shown that these thermal, statistical states also include a fairly wide set of decoherence free, multi-DFS states. 
Under this extended perspective, conditions (ii) and (iii) become particular instances of more general requirements 
for decoherence free propagation, and their discussion is covered under the general case. 



C. The finite temperature coherent-product states 

The density matrix generalization of the Davydov coherent-product ansatz ]27| that we purport to employ, can be 
regarded in the current context as arising from a particular set of initial density matrices /5(0). Following this point of 
view, we introduce the coherent-product distributions in a self-contained manner [independent of the approach in ref. 
p7[], via a formal protocol for the preparation of the initial state. An alternative possibility specific to the spin-boson 
model is also suggested in Sec. VI. 



1. Preparation of the initial coherent-product state 



Consider thus the following procedure: 

i) Let the register couplings (^) be modified by external means such that the original structure of the register DFS 
is refined, while the strength of the coupling eigenvalues is arbitrarily altered. For notational convenience, we assumc 
that the DFS degeneracy is completely lifted, but more realistic cases can be easily substituted by allowing some of 
the coupling eigenvalues to coalesce at identical vàlues. In other words, if is an orthonormal register basis 

yielding an irrep on 7Z for t > 0, let all DFS be reduced to the 1-dimensional subspaces of the individual states Uj> R v ), 
with corresponding projectors Pj u = {<(>%) {(f'jvl· ^~fik,Ju denote the new coupling eigenvalues, decomposition (|) for 
the register couplings becomes 



J,v 



and the corresponding Hamiltonian amounts to 



H' = Hn 



5>- 



Jv 



E ^k,Ju (Xkqbl + Xkqbq 
k . q 



ii) At some time to < the register, in an arbitrary state p R , is brought into contact with the uncorrelated 
environment in thermal equilibrium at temperature T, such that 



p(t ) = p 



R - n B 

Pt 



(19) 



with 



Pt 



Z b 



exp 



kfíT 



iii) The total system is subsequently allowed to evolve under the modified Hamiltonian H' above, until time t = 
when the original couplings are instantly restored. Using the appropriate version of Eq.(|^), the state so prepared at 
t = is easily seen to be 



,5(0)= J2 Pj»Vj»P R Pb,t{V?,„,) ] Pj> v > 



(20) 



Ji/.J'u' 
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where irrelevant phase factors have been absorbed into the register basis states |</>j„) , and we introduced the short-hand 
notation 



yf v = ex P 



exp 



i 

h 



XI huJ q b\b q + ^ (pju,qt>\ + ÏÏju,q b q 



1 1 

Also, the adimensional parameters (5 Ju _ in the last line above are defined by 



(21) 



Xkq 



The environmcnt contributions 



Vpvf>B,T {yp'v')^ to tne density matrix ( p0| ) can be brought to a more convenient 



form if the factors Vf v are written as 



% - Huj q b % É ° 



exp 



i2_^uj q t \!3j Vfí 



exp 



^ Pjv,q b q — Pj V! qbc 



1 



• exp 



■ exp 



exp 



The last three factors on the right hand side can be easily rearranged as exp 
to yield 



l - twqblb^j t a 



Vf v • exp 



X hujgblbq t 



exp [i9j v ] ■ exp 



£(/W°) 6 í-/W Q ) 6 «) 



(22) 



where the displacement parameters f3j v , q (0) = 0j uq (é lu ' qto — l) can be retained as final parametrization variables. 
Use of identity ( |22| ) in expression (^Ü|) now leads to the following defining expression for the sought class of initial 
density matrices: 



p(o)= p ^ ex p 



Jv,.Tv> 



£ 0WO)&J - #W )M 



L 9 



exp 



p 



JV 



(23) 



where the phase factors 9 j v are absorbed again into the corresponding register states. 

The density matrix ( |23| ) describes a register entangled with an environment in a superposition of Gaussian states, 
as shown by the environment reduced density matrix 



p B (0)=Tr R [p(0)}=j: TrRÍ ^ R) «1» 



L·j2hw q (bl-0* JV (o))Q>q-0j»(o)) 



k R T 



In the particular case when all displacements /3ji,(0) vanish, or alternatively, when the preparation procedure is 
skipped (to — 0), the environment remains in a thermal, uncorrelated state. Another interesting limit anses when the 
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environment is cooled at zero temperature, such that p® = |0b)(0b|, and the register state p R corresponds to a pure 
state, p R = |<Í>· R )(<Í> iï |. Then the total state (|2^) also describes a pure state, of the form 



where the coeficients are given by cj„ 



l*> = E c ^l^)®K^(o)}) , 



Ju 

b%\$ R ) 



Tli\ c Jv\ 2 = li an d the environment factors read |{/?ji/(0)}) 



exp 



|0fl). One recognizes without difficulty a coherent-product Davydov state. Inter- 



estingly enough, the general expression (23) can be brought to a form reminiscent of the Davydov ansatz as well, by 
the formal artifice of square-root factorization. Indeed, if we introduce a square-root representation for the density 
matrices p R and p R , i.e., 



1 K ■ (7 K ) 



the entire density is factorized as 



with 



P(0) = 7(0) • t+(0) , 



7(0) = £751(0)^731(0) 

Ju 



The register and environment factors in the decomposition above now read respectively 
and 



7|„(0) = exp 



^2(0Ju, q (O)bl-j3*j Víg (Q)b q 



•7t 



The total density matrix p(0) acquires in this way a structure similar to the density matrix for a Davydov pure 
state, where the various pure state factors are replaced by [nonhermitian] operator factors. Moreover, the register 
factors are orthonormal in the sense of the usual trace inner product on the space of linear operators on 1Z, i.e., 
(7jt(0)l7jv(0)) = Trii[jj u (0) ■ { A fjt v i{Q)y] — àjv,J'v<, and the environment factors 7^(0) are easily verified to be 



square roots of Gaussian distributions of the form exp 



This formal analogy 



prompts the designation of states described by Eq.(|23j) as density matrix coherent-product states. It was found in 
ref. p7[ , and it is shown shortly for the case at hand, that such states preserve their form throughout an evolution 
driven by a Hamiltonian of type provided the register states are restricted to a direct sum of DFS. 

Before proceeding in this direction, let us stress that the register density matrix p R in expression ( [23] ) above is not 
identical to the initial reduced density matrix p R (0) unless all displacements (3j v . q (Q) vanish, and the total density 
matrix reduces to the thermal, uncorrelated form (^9|). In all other cases, the correct reduced density matrix at t = 
is entangled with the environment, and reads 

p R (0)= £ r,J»,J>v>(0)PjuP R Pj>»> , 



where 



exp 



0Wo)&J - /W°) 6 * 



Pb,t exp 
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Thus the matrix p R represents only a register parametrization variable for /3(0), just as the displacements f3j Uiq (0) 
parametrize the environment contribution. More precisely, the correlation amplitudes between the register states 
\<j)%) projected by Pj„ are given as 



and include both an environment-mediated component contributed by r\j v ,j' v i (0) and a parametric register component 
represented by (4>%\p R \4'%) ■ 



2. The time- dependent coherent-product solution 



The density matrix evolved under Hamiltonian (||) from initial density matrix ( p3| ) can be calculated straightfor- 
wardly with expression (||) for the evolution operator. After a rearrangement of the type leading to identity (]2^), it 
is found to read 



p(t) = cxp 



Ju : J' ' v' 



exp 



~H R t 
n 



(24) 



where the phase factors Ojv(£) can be shown to amount to 



(25a) 



and Uj v (t) are environment displacement operators of the form 



Uf„(t) = cxp 



(25b) 



In Eqs.(|25|) above, the environment displacements 0j v ,q(t) evolve as 



,„,,,„: , íW0) + £ii) e -..-g 



(26) 



As anticipated, solution (|2J) indeed preserves the coherent-product character of the initial condition (|23|). In a 
square-root decomposition p(t) — j(t) ■ J^(t), a representative [nonhermitian] square-root j(t) takes the form 



with factors given by 



Ju 



exp 



-H R t 



7^(0) , 



and 



WuWW'At)) = Tr R [*jtf v (t) ■ (j$„,(t)V] = 8 Ju ,j, v , , 



íjÀt) = exp 



■7l 



Obviously, the environment remains in a nonstationary statistical superposition of Gaussian states, with a reduced 
statc 
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\t)=Tr R [p(t)]=Yl 



Tr R {P JvP R ) 



exp 



Jv 



It is worth noting in particular that Eq.(p4[) ofFers a closed density matrix form for the well-known exact solution to 
the decoherence problem with an uncorrelated, thermal initial state of type ( |Ï9|) |3C| l. From the above observation 
on the state of the environment, we obtain as a direct corollary that the evolution of an initial state under 
Hamiltonian (Qj drives the environment into a nonequilibrium state of type $2$), strongly entangled with the register, 
unless the initial register state is restricted to a single DFS. Conversely, and with reference also to the preparation 
procedure described earlier, the coherent-product solutions can be regarded as straightforward generalizations of 
the Standard, thermal solution. The formal essence of this generalization is seen in expression ( p6| ) for the time- 
dependence of the environment displacements Pj v<q (t). That is, the coherent product form ( |24| ) for a generalized 
initial state (|2^) differs from the coherent product form for the thermal solution solely by the presence of nonzero 
initial displacements /3ji/, g (0). It can be said therefore that the Standard solution is generalized here by allowing for 
nonzero initial environment displacements in its coherent-product expression. 

Let us point out also that in the zero-temperature, pure state limit, corresponding as before to p R — \§ R )(<fr R \ and 
p R =0 = |0s)(0b|, the coherent-product density matrix ( |24| ) describes a Davydov-type pure state evolving as 



!*(*)) = exp 



-H R t 



E 



cj v e 



Wu)®\M)) > 



with cj v 



exp 



E{0Ju, q (tH-fiv, q (t)b q ) 



2 = 1, and the environment coherent states \f3.j u (t)) defined by |/3j„(í)) = 
|0s). As must be expected, we recover the Davydov state example ( |Ï7j ) of Sec.IIIB. 



The apparently different phase factors are due to the representation of the environment coherent states in terms of 
unperturbed modes b q , rather than displaced, DFS-specific modes bj q . 



In the next Section we show that the class of thermal coherent-product states can generate multi-DFS decoherence 
free reduced register distributions, under reasonably relaxed constraints on the coupling constants pj q and on the 
symmetry of the environment. The general form of the reduced register state reads now [compare to Eq.(pl)] 



p R (t) = Tr B [p(t)} =exp 



-h HRt 



Jv.J' v' 



■ exp 



n HRt 



(27) 



and shows that the effect of the environment is conveniently concentrated in the bath-mediated correlation factors 



Wv , (í) = e m.,^,,m TrB \uf u (t)p BíT (uf,At)) 



(28) 



where the phase factors Qj v and the unitary displacement transformations Uf v are defined by Eqs.(P5[). The trace 
in Eq. (^8|) can be readily calculated, either in the usual manner, via the symmetric order generating functional for 
the harmònic oscillator or by the techniques of thermofield dynamics P7| . It is also helpful to use the prior 
observation that expression ( |2§| ) differs from the corresponding result for uncorrelated, thermal initial conditions only 
by the specific time-dependence of the displacements f3jv,q- The ensuing expression for the bath correlation r]j v ji v i 
is of the form 



where the phase factor $ jvj'v' amounts to 



*Jv,J>W = [PJVA (*) ft>v> t9 (*) ^ (*) PJ>*« (*)] 



(29) 



(30a) 



while the dissipative factor r 1 j v ,j'v' is given by 



r w (í) = \ E \Pjv, q (t) - Pj>u,, q (t)\ 2 coth (^Ç) 

q ^ ' 



(30b) 
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For the spin-boson model ([j]) with an uncorrelated, thermal initial condition, it can be verified starightforwardly that 
these general forms reduce to the expressions recently derived in rcf . [p9| . 

The physical becomes apparent when the trace in Eq.(p8|) is expressed as 

a thermal average, to the result that 

r? W (t) = JPMO-esAO) J2 eXH ~ Eí ^ }/kBT \ {n q Mt)}\{n^At)}) , (31) 



where 



\{n g ,ju(t)}) = exp 



IK» , 



with |{n 9 }) the excited bath state with n q quanta in mode q, and E{ n y = Y] Tiuj q n q the excitation energy of the state 

i 

\{n q }). Recali that the states \{n q .j^(t)}) can be understood in terms of the displaced modes 



52(f3j v , q (t)bl-í3*j v , q (t)b q 



■bl- 



exp 



bl-ftvM) 



as [timc-dcpcndcnt] environment states with n q displaced quanta of mode q excited over a displaced vacuum 



|{/?,,j„(t)}> = exp 



10} 



That is, 



\{n q ,Mt)}) = II^Ï (í'UW)"'' 



Hence the bath factors r\jv,j'v' (t) represent, up to a phase factor, the thermally weighted sum of interference ampli- 
tudes between the time-dependent coherent environment states entangled with the register states |<^j„) and \(f>j/ v /) 
[i.e., with the register factor Pj v f>R • Pj'v 1 ] in the overall dcnsity matrix. 



IV. DECOHERENCE FREE STATES ENTANGLED WITH A NONEQUILIBRIUM ENVIRONMENT 
GENERAL CONDITIONS IN THE COHERENT-PRODUCT ANSATZ 



A first hint that nontrivial, multi-DFS decoherence free states are possible in the thermal coherent-product ansatz, 
comes from the straightforward observation that the harmònic time-dependence law for the environment displacements 
given by Eq.(E6|) admits the stationary points 



7S - fiJq 



If the initial state is such that the environment modes are displaced exactly over these positions, the reduced environ- 
ment state remains stationary and the bath-correlation factors in expression (p7| ) for the register reduced state varies 
only through the phase factors Qj u . The origin of this effect can be seen by noting that in this case the action of the 



displacement operators (|25|b) on the excited states of the unperturbed environment generates, up to a phase 

factor, the similarly excited states of the corresponding DFS-specific Hamiltonian, i.e., 



Uf v \{n q }) 



exp 



E 



hui„ 



\{n q }) = \{nj, q }) 



where [see Eqs.(|Ï0|)] 
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HBj\{n,j. q }) = 



V 8 



Taking this into the environment factors 



of the coherent-product density matrix (124) , shows 



U? v (t)PB,T(U?, u ,(t))[ 

that the register states are actually correlated with eigenstates of the environment Hamiltonians Hb,j, which vary 
only through phase factors under the action of the total propagator (flï|). 

Returning now to the bath-correlation factors in the register reduced state (p7|) , observe that the remaining phase 
factors Oj„ become proportional to the energy-shifts (|ïo|c) and linear in time, as 



Consequently, the evolution of the register state (|27]) acquires the unitary, quasi-unperturbed form 



p R {t) = exp 



j-(n° n° )t- p ~R. p 



exp 



-Hnt 



exp 



£ Vj v ,j'u'PjuP R Pj'w 

Jv,J' ' v' 



• exp 



As in the first pure state example of Sec.IIIB, an exact and nontrivial unperturbed propagation requires that the 
reference energies of the contributing DFSs be identical, such that 



flLü q ^— ' ?ÍU! q 



= , 



and that some correlation factors between distinct DFS [J ^ J'\ be nonzero. Otherwise, as will be shown to happen 
for model ([!]) in an environment with a one-dimensional spectral density, the correlations reduce to t}j v ^ji v i = 
and the presence of the bath is effectively erased in p R (t) . The latter becomes a trivial block diagonal distribution 
on disjoint DFS, and individual decoherence-free states become eventually pointer states p2| . 

Let us now seek the general conditions for the unitary and unperturbed propagation of the reduced register state in 
the coherent-product ansatz, starting from the requirement of time-independent bath factors, i.e., 



dt 



fdQ Ju dQ 
\~~dt 



J'v> 



Jv,J'v' 



dt 



dt 



dT 



Jv,J>v' 

dt 







(32) 



Substitution of the corresponding expressions, Eqs.([25|a) and fl3C|), including the explicit time-dependence ( |26|) for the 
bath displacements, and a little àlgebra yields a Fourier-like sum of the form 



d 1 \ -> / 

dt lli7lJ ^ J ' u ' = F ° + 2 ^ ' í «' +í 



(33) 



with coefhcients 



^ = £ 



huj n 



hu> n 



(34a) 



F q ,+ = (»*j q 0%, q ~ ti> q Pj> v >, q ) + (iï q 0j> v ', q - ^A,) - coth (Jj£g) ~ 0j> v >, q ) " fi> 9 ) , (34b) 

and, respectively, 
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F q ,- = (ti q P°ju, q ~ ti' q P%, q ) + " Vj'A,) + coth {J^ ) (p\ q - (3% v , tq ) (n* Jq - . (34c) 



r, 1 í , d\q\ 

F + - du- " 




Here the labels J'u' and Jv have been dropped for notational economy, and we have denoted 

The Fourier-like nature of the right hand side in Eq. (|3^) can be easily exploited under the customary assumption of 
an d-dimensional environment with a continuum spectrum, oj q — > cu(q). We also assume, without loss of generality, a 
spherical dispersion relation of the form u(q) = o;(|g|), and an environment volume set to unity. In this case, the sum 
in Eq.(|33"|) can be rearranged into a Fourier integral by subsuming contributions from degenerate modes on spheres 
S Uq =ío, to the result that condition (|3^) is brought to the form 

2 j - düj , T dS q F* j+ | +(<f dS q F q ,-\ = . (35) 

Since the cancellation of a constant term Fq ^ demands 

dS q F q .±J ~ í(«) 

and, in all likelihood, divergent displacements for the static ground mode w = 0, it is seen that an unperturbed 
evolution of the register requires Fq = 0, i.e., 

dq \H±t = í dq \t^£ , (36a) 

ÏWJq J tlUJq 

as well as 

<j)dS q F q: ±=0 (36b) 

for all Jv and J'v' . 

Condition (|36|a) is easily recognized as the energy-shift condition (i) [Eq.(|Ï5|)a1 for the pure coherent-product states 
of Sec.IIIB, and also retrieved for the particular example in the beginning of this Section. Conditions (|36|p) reduce 
eventually to the simpler linear system 

dSA, g (iï q - fi> q ) = , VJ, J', v , (37) 

when noted that according to the explicit expressions (Q) the terms in coth (huj q /2kBT) contributed by the dissipative 
exponent Tj v ji u i cancel separately. Since this system is temperature independent, it necessarily applies also to the 
pure state case of Sec.IIIB. In this limit it is seen to ensure that the evolution of the environment coherent states 
preserves their overlap. By analogy, in the finite temperature situation it secures the synchronization of the Gaussian 
environment distributions, so that their thermally averaged interference yields a time-independent bath-correlation 
factor rjj Ut j' U > . 

Thus the unperturbed evolution of a register distribution in the generalized coherent-product ansatz requires si- 
multaneously: 

i) identical energy reference points ( energy-shifts ) for the decoupled environments generated by contributing register 
DFS [Kqs.^fi) or 

ii) a specific synchronization (phasing) of the entangled environmental modes, through proper displacements fij v , q 



[Eqs.(3VJ; 



For any register of type ([!]) states compatible with conditions (i)-(ü) always exist, because for any DFS characterized 
by \ij q there exists another DFS characterized by (—(J,j q ), and related to the former by a reversal of all qubits along 
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direction z [| T)« ^ | |)nj I l)n - * I T)n], in every member state. This ensures that condition (i) can be satisfied. 
Condition (ii) is satisfied at least by the trivial solution to Eqs.(|37|), (3j vq = 0ji v i „ = f° r au li which corresponds 
to the case of stationary displacements discussed in the beginning of this section. Nontrivial solutions are essentially 
conditioned by the degeneracy of the environmental modes. It is worth noting that the number of distinct DFS 
contributing to a decoherence-free distribution of the type discussed here is theoretically arbitrary [if finite], unless 
the environment is 1-dimensional and displays only two-fold degeneratc modes [uj~ q = uj q \. For the lattcr case, the 
number of distinct DFS involved cannot exceed 2, since the number of unknown j3j v . q -s in system ( |3^ ) must exceed 
the number of constraints. 

However, in addition to conditions (i)-(ü) above, the existence of nontrivial environment-entangled state is essentially 
limited by the requirement that [some of] the bath correlations r f]j v ^,j i v l must be nonvanishing, which means that the 
corresponding dissipative factors Tj u ji u i must be finite, Tj v> ji v i < oo. Assuming a nontrivial solution to system 
( |37| ) does exist, substitution of the resulting bath displacements in expression (|30|b) yields straightforwardly [the 
time-dependent terms vanish by conditions (|36ï) 



íj v ,j>v> oc- düj— — coth 

2 / OU) 



htv 
2k B T 



dS a 



+ 



(J-Jq ~ HJ'q 



tluJ 



(38) 



It follows that Tj Vt j> v i cannot be finite unless its value Tj j, for stationary initial displacements \0j vq = /3j<„< „ = 0] 
is also finite, 



r °/,./' « \ J dw ~ d ~^ huj ^ 2 coth C^t) / dSq ^ Jq ~ MJ ' ? I 2 < 00 



(39) 



This intrinsic stationary dissipative factor r j j, is seen to be characteristic of the two contributing DFS, and not 
of individual contributing states. It sets an upper limit on the magnitude of the corresponding bath correlation 
n)Jv,j'v'i an d therefore on the amplitude of correlation between the associated register states [see Eq.(|^) for the 



reduced register state p (t)]. It also increases with the temperature 



oo as T — > oo] regardless of the exact 



density of states or of the form of the coupling constants iij q , and gradually shrinks the set of register distributions 
compatible with unperturbed propagation toward trivial states, block diagonal on DFS [rjj^j'v' — > as T — > oo for 
J 7^ «/']. Moreover, the presence of properly phased coherent oscillations of the bath modes „ — Pj v , q \ > 0] 

results invariably in decreased bath correlations r)j v .ji v i and decreased cross-DFS matrix elements in p R (t). Hence 
the preservation of an unperturbed register evolution seems to involve a trade-off between the stabilizing action of a 
nonequilibrium environment and the magnitude of the conserved register correlations. 

On the other hand, the finite or infinitc character of T j Jr at finite temperatures does depend on both the density 
of degenerate bath modes, and the coupling constants for the input DFS. To find the source of this effect, let us 
examine the thermal average representation ( |3l"| ) of the bath correlation r\j v j' v i in the limit situation when Eqs.(|3^) 
are trivialy satisfied for 0j, v , A = f3% q = [Pj v , q = -{n. Jq /huj q ), 0j> u >, q = -(jijt q /hu) q )], T Jv j, v , -> Tjji = T 7 j, and 
f]Jv,j'u' -> VJJ' ~ ex P[- r j,j'í- To this end, substitute (3 Jv , q -> f3j tQ = -{pL Jq /Tw q ) and (3j> v > tq -> f3j> >q = ~{nj' q /hu> q ) 
and write the resulting overlaps ({n qt j}\{n q ,j'}) in the form 



(Kj}IK^'}> = ({«Jlexp 



E 



bl 



hiüq q 



huj„ 



exp 



Q 



hiÚn 



■bl 



híú a 



IK» 



ò^exp 



Vjq Mj'q - 

huj„ 



exp 



(40) 



where ctj j> = exp 



(V 2 ) E {VjqVJ'q ~ HjqH*J> q ) I {^qf 



in the latter expression can be calculated as 

Qn Q 



exp 



:i/2)eimj ç -mj' 9 i 2 /(^ 



The vacuum averages 



(-1)' 



dX n d(X*) 

with P n a polynomial expression, to the result that 



<0, |exp [-X*b q ] exp [Xb\] \ q ) = P n (A, A*) exp [~|A| 2 ] 
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({ n q,j}\{ n q,J'}) OC (0j|0j/) = CXp 



and, equivalently, 



exp [-T° Jt AT)] cx (0j|0j,) = exp [-T% jr (T = 0)] 



(41a) 



(41b) 



As in Sec.IIIB, here |0j) denotes the displaced environment vacuum created by the J-th DFS. 

Similarly to the pure state case therein, the finite character of the stationary dissipative factor T°j j,(T) is 
seen to be conditioned by the unitarily equivalent or inequivalent character of the environment vacua gener- 
ated by the register DFS, or equivalently, by its zero temperature value T°j j,(T — 0) [this latter form of the 
condition can be obtained directly from expression ( |39| ) for T j Jt (T) in the zero temperature limit] . Clearly, 
the displaced vacua are unitarily inequivalent [have vanishing overlap], and T°j j,(T = 0) also vanishes, when 
(0q,j\0 q ,j>) oc exp[-(l/2)|jtij g - n J/q \ 2 /(huj q ) 2 ] — > as q -> [\(jij q - fij> q )/huj q \ diverges as g -> 0], unless the 
density of modes compensates for the contributions from low-frequency modes. At finite temperatures the effect is 
furthcr amplified by the thermal excitation of low-frequency states. On the other hand, the situation can improve 
considerably when \{fJ,j q — fij' q )/huj q \ remains finite as q — > 0, and the displaced vacua remain in unitarily equivalent 
Fock spaces. 



We may conclude that decoherence free states involving múltiple entangled DFSs are possible in the thermal 
coherent-product ansatz if conditions (i)-(ii) above are simultaneously satisfied alongside condition 

iii) the stationary part jr of the dissipative factors Vj v ji v i must be finite; a necessary prequisite is that the envi- 
ronment displaced vacua associated with the contributing DFSs belong to unitarily equivalent Fock spaces. 

Conditions (i) and (iii) impose limiting restrictions on the register-environment coupling, by demanding both a 
particular structure of the coupling constants, and a suitable density of environment states. Condition (ii) can be 
always satisfied for continuously degenerate environment modes, by a proper choice of initial displaccments. Finite 
degenerate modes [e.g., the 2-fold degenerate modes of a 1-dimensional environment] can also comply with condition 
(ii), provided the number of distinct displacements is appropriately limited. 



V. EXAMPLES: THE SPIN-BOSON MODEL 



Let us illustrate the arguments of Sec. IV on the specific case of model ([j]). Under a common functional prescription 
for the density of states, it turns out that a single qubit register [N = 1, \qn = x( w <?)] can display a finite stationary 
dissipative factor T° only in a 3-dimensional environment. Not surprisingly, a similar behavior is also seen in multiqubit 
registers with weak collective decoherence [N > 1, Xqn — xi^q)]- However, linear registers with individual decoherence 
[Xqn = xi^q) exp(iç • r n )] present states for which exp[— (l/2)\/j,j q — f-j'ql 2 / ( nLü q ) 2 } remains finite as q — * 0, and which 
generate finite r°-s even in a 1-dimensional environment. 



A. Single-qubit register 

For a single qubit register [N = n = 1] let the model coupling constants assume the form 

Xq = X{Uq) ■ 

There are only two one-dimensional DFS, such that Jv =| or Jv =J., and the associated couplings jij q become 



/if q — —(Jti q — xi^q)- Assuming a total coherent-product state, the reduced register state (27) has the form 

pK(t) = e"*-* 4 • [| t)(T \p r \ TXT I + I IXI \p r \ IXI |+ 
+v(t)\ TXT \p r \ IXI I + »7*(*)| IXT \p h \ TXT \p r \ l}*(T I] ■ e^ É , (42) 

with the environment driven correlation r/(t) is given by Eq.(p^) with appropriate substitutions. 
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The state p R {t) becomes decoherence free provided rj(t) = const. under the corresponding conditions (i)-(iii) of 
Sec.IV. Condition (i) [Eq. " 
(ii) [Eqs.@] reduces to 



6a)] is satisfied by default, due to the particular form of the coupling constants. Condition 



ds q p° u = i ds q p° Lq = 



(43a) 



whilc condition (iii) for a finite stationary part of the dissipative factor T requires [see Eq. ( |39|) 



dS q \2 X (u)\< 



í d JM G{uj) 

J diu 



coth 



2k B T 



< oo , 



(43b) 



,d-l 



where G(u>) is the density of modes at frequency lo. 

The density G(u>) grows as area(S u ) — § dS q uj , 

S u 

characterized by a natural ultraviolet cut-off frequency u> c , 



where d is the dimension of the environmcnt, and is 

which sets the upper limit for the rate of dissipation 

d-l . 



processes in the environment. This necessary feature is usually accounted for by setting G{u>) oc uj d ~ Y exp[— uj/u c ]. It 
is also common to assume a quasi-linear dispersion d\q\/duj ~ const. and a coupling x(lü) oc \fw. Alternatively, one 
can transfer the cut-off exponential onto the interaction strength, and set |x(w)| 2 oc wexp[— u>/u) c ], while leaving the 



density of states as G(u>) 



uj d 1 . In either case the final prescription reads 

(d\q\/diü)G{u;)\x(Lü)\ 2 = \(h 2 u c /2)(u/uj c ) d exp[-ij/ij c 



where the scaling factors are chosen such that A is an adimensional constant. In the following we favor the second 
interpretation, which can be straightforwardly extended to cases with couplings unisotropic in the wavevector q. 



For the present single-qubit case, \\(uj)/huj\ 



oc to 



exp[— lü/lü c ] — > oo as q — > 0, and the stationary dissipative 



factor r° cannot be finite [equivalently, the displaced environment vacua associated with each of the two qubit states 
cannot be unitarily equivalent] unless this singularity is balanced by the density of states. This is also evident from 
the integrand in Eq. (43 d), which behaves in the low- frequency limit as uj d ~ 3 . In particular, a 1-dimensional density 
of states proves insufficient to counteract the low- frequency contribution even at zero temperature, i.e., the displaced 
environment vacua are unitarily inequivalent. A 2-dimensional density suffices at zero temperature, implying that 
the displaced environment vacua are unitarily equivalent, but fails at finite temperatures. Only in a 3-dimensional 
environment does T° acquire finite vàlues at both zero and finite temperatures, meaning that the displaced vacua 
are unitarily equivalent and the generated DFS correlations are stable under thermal excitation. In this case, the 
exact expression of T° can be given an analytical form in terms of the generalized Riemann zeta function £(<?, z) — 



(l/f(g))Me- 





where r°(0) = A 
Expression 



/(l 



; ) [here f (g) denotes the Gamma function] as 



r° 



TlLür 



= r°(o) 



k B T 

TlLür 



C 2 



k B T 

TlÜJr 



- 1 



(44) 



shows that the stationary factor T° increases monotonously with the temperature, behaving as 

r° w r°(Q) 1 + (tt 2 /3) {k B T/hu c f in the low temperature limit [k B T « huj c ] and as T w 2r°(0) (k B T/Tiuj c ) at 

high temperatures [k B T >> %u) c ]. Accordingly, the maximum amplitude of register correlations (rj ■ (| \p R \ J,)) that 
can be carried in an environment-entangled unperturbed propagation decreases fast to zero in a classical environment, 
with temperatures k B T > hoj c . A nontrivial unperturbed propagation is seen to require a low-temperature, quantum 
environment. The widest accessible range of register correlations is supplied therefore by zero-temperature, Davydov- 
like pure states of the joint register-environment system. For environment displacements /3|(j.) g (í) satisfying the 
phasing condition (f43|a), it can be checked that such states acquire the form 



|*(t)> = e-* 



-íll 



c t |T)®I/3t(í))+c;U)®I^(í))] 



wi t h I^TU)(*)) = ex P 



E(/3m) q (t)b\-(3 m)q (trb q ) 



\0 B ), and (/3t(*)|^(í)) = </3 T (0) |/3 a (0)) . It was also taken 



into account that under condition (43i) the phase factors 8|(n(í) reduce to 9|(í) = Oj^(í) = 4ííí, with 
fi = - Jduj(d\q\/cL·)G(uj)[\x{uj)\ 2 /hLj}. 
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B. N-qubit register with weak collective decoherence 



The situation is vcry similar in the closely related case of a multi-qubit register [N > 1, < n < (N — 1)] with weak 
collective decoherence, when again 

Xqn = Xi^q) ■ 

The register DFSs are now indexed by the total spin projection on the z-axis, e.g., 

JV-l 

j = E ^ . 

where si"' is the eigenvalue of labcling the z^-th state of the J-th DFS Sj, and the corresponding coupling 
constants read 



Due to condition (|36|a), the register states contributing to a given environment-entangled decoherence- free distribution 
can only belong to either of a pair of DFS, {Sj U<S_ j}, characterized by coupling constants J • x(w 9 ) and (— J) • xi^q), 
respectively. As a result, an entangled decoherence free state is of the form 



p R (t) = exp 



+VJ^-Jv>(t)Pj„p R P~jv> +V*ju,-j V '{t)P-Jup R Pjw ] - exp 



n HRt 



(45) 



where the bath correlations between states in the same DFS are explicitly shown to be time-independent. The cross- 
DFS bath correlations rjj u .~j V '{t) become time independent provided the environment displacements are such that 
[see Eqs.©] 



fdS q fi Vjq = fdS q lfi_j Vtq = Q,Vi> 



(46) 



and remain nonvanishing if the stationary dissipative factors rj_j are hnite, which amounts to 



diu 



■ coth 



/ Tiuj 



\2k B T 



dS q |2 Jx(u)\ 2 = J 2 T a <oo . 



(47) 



In Eq.(^) above, T° is the stationary dissipative factor discussed previously for a single qubit. Evidently, under 
the same prescription for the density of environmental modes, a similar discussion applies and shows that the bath 
correlations r\j v ,~j V ' can be constant and nonzero only in a 3-dimensional environment. Note, however, that thcir 
magnitude decreases strongly for register states with spin projections |J| > 1. In the zero-temperature, pure state 
limit, decoherence free register distributions are generated by Davydov-like states of the total register-environmcnt 
system with at most dim(Sj) + dim(S-j) — 2 x dim{Sj) distinct terms. 



C. N-qubit register with individual decoherence 

Consider now a N-qubit register [N > 1,0 < n < (N — 1)] in an individual decoherence regime, characterized by 
coupling parameters 

Xqn = X(^q)cxp(iq ■ T n ) , 

where r„ denotes the position vector of the n-th qubit. We also assume, for convenience, a linear geometry, such that 
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The DFS of this register are trivially 1-dimensional and correspond to the eigenstates |{sj^}) of the unperturbed 
Hr, such that the coupling constants jij q now amount to 



JV-l 



,inqaM 



n=0 



Note that the DFS index J is no longer associated to the total spin projection as in the previous case, but labels 
individual unperturbed eigenstates. 

Select a specific eigenstate |{sj™^})- In order to detect a set of entangled decoherence free states involving |{s^}), 
consider first the corresponding condition (j36[a), which limits the set of compatible coupling constants fij q . A straight- 
forward solution is to seek states for which \nj q \ = \^j aq \,^q- Given the linear translational symmetry of the register, 

it is not hard to observe that the set of such states includes the N states \{sj }) [0 > m > (N — 1)] rclated to |{sj* }) 
by a cyclic permutation modulo N, i.e. 



„(») _ .(»-««) 



Vn 



Indeed, the corresponding coupling constants differ from fij oq only by a phase factor, respectively 

iqma 



AU, 



a'l 



A coherent-product total state involving the N register states thus selected produces a reduced register state [see 
again Eq. (EJ 



p R (t) = exp 



-H R t 



VJ m ,J m >{t)Pj m P R Pj ml - exp 



-H R t 



(48) 



for which the environment mediated correlations ï?j m ,j , (t) become constant in time when the bath displacements 
comply with [Eq.@] 



(49) 



Nonvanishing f]j mi j ,-s require in addition a finite corresponding dissipative factor, i.e. [Eq.(B9)] 

2 



r o J oci [ düJ d M coih (^_ 
JmJ ™' 2J du \2k B T 



dS a 



< oo 



(50) 



Let us assume again the usual prescription for the density states, as described for the single-qubit case. If we also 
introduce the trànsit time t s , defined by qa — u)t s , it can be seen that for any m and ml the integrand under the 
surface integral above reads 



N-\ 



\(l· l J m q ~ l· l J m >q)/ huj q\ 2 = M(x(^)/^) sm[(m' - m)u)t s /2] ^ 



ots „(")|2 



and behaves as u> when q — > [recali that |x(w)| 2 ~ lo exp[— lü/lu c ] ]. Hcnce the entire integrand in the stationary 
dissipative factor r° j l behaves in the low-frequency limit as area(S u ) ~ uj^ 1 [for d\q\/du> ~ const.] and j 
remains finite in any environment, at any finite temperature. Equivalently, the displaced environment vacua cor- 
responding to the register states \{sj }) and [{sj™" 1 ,}) are unitarily equivalent and the corresponding correlations 
between the register states are thermally stable. 

One can also consider the mirror permutation states with 

„(«) _ AN-n) _ (N-n+m) 



and 
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i(N —m)qa 



PJgq 

xK) 



which are compatible with condition (36 a) as wcll. The cquilibrium dissipative factor of any pair of such states reads 
r^- - = r ° Jm j / , while the factor corresponding to a direct permutation state and a mirror permutation state is 

given by 



ii / , d\q\ , / fko 

r T -j oc 2 / diu-P- coth 1 

•>mJ„,> I du) 



\2k B T) J dSq 



i \ N - 1 
X M \- (n) . 
— — > s\ Sil 

n=0 



m + m! — N\ 
n H 1 ut s 



(51) 



and again proves to be finite in any environment. It follows that any environment-entangled distributions which 
involve the 2N direct and mirror permutation counterparts of |{sj™ }) and satisfy conditions ( |49| ) generate nontrivial 
register mixtures propagating in an unperturbed manner. Note however that in a 1-dimensional environment such 
a mixture can only accommodate two distinct pure states [see Eqs. (ft9[)]. It is also worth noting that in this case 
the precise magnitude of the stationary dissipative factor T° is strongly dependent on the input states. The simplest 
example of a joint register-environment state yielding a decoherence free register evolution is again a Davydov-like 
pure state, which in the present case takes the form 



!*(*)> = cx p 



--Y 



eaí n H 



-m | « j 



® \Pj m {t)) +C r 



i4 n) )®i/?j m (í)) 



The expression above is obtained for coherent environment states \(3j m {t)) and \ j3-j (t)) with displacements (3j m ,q{t), 
j3j q (t) satisfying synchronization conditions similar to Eq. (ff9|) . In this casc the overlaps (f3j m (t)\f3j m , (í)), 
(f3j (í)|/3j ,(í)) and {f3j m {t)\(3-j r (t)) are constant in time and the phase factors Qj v (t) [see Eq.(^a)] become 
Qj m (t) = 0- (t) = 6j (í) = (l/h)Qj tym, with Slj a constant energy shift appropriately defined. 



VI. ELEMENTARY ENTANGLEMENT REPHASING THROUGH QUBIT FLIPPING: 
CONTINUOS-TIME COHERENCE EVOLUTION DURING BANG-BANG CONTROL 

The coherent-product ansatz underlying our discussion of multi-DFS decoherence free states was introduced in 
Sec. III via a formal preparation protocol relying on an explicit alteration of the register interaction with the envi- 
ronment. But a practical approach to the preparation problem may not necessarily require direct manipulation of 
the environment. As a counterexample, let us point out that at least one method of decoherence control available 
in the literature involves an elementary rephasing of entanglement of the type discussed here. We refer specifically 
to quantum bang-bang control inspired by the multi-pulse decoupling techniques of NMR, which counteracts 
decoherence through a train of identical spin-flip cycles. Each cycle generates a revival of coherence through a pair 
of coherent 7r-pulses that alternately flip the state of the register. When this process is examincd in the context of 
the extended density matrix solution of Sec. III, the cause of this revival effect is distinctly exposed as an elementary 
adjustment of the entanglement with the environment. In other words, spin flipping provides a simple working pro- 
cedure for manipulating entanglement. In a supplementary outeome, the exact density matrix solution brings forth 
an alternative picture on the source of the refining power of so-called symmetric cycles ||Ï2| . The latter introduce a 
straightforward adjustment of the bang-bang technique, which lowers the working cycle frequeney by a factor of at 
least 2. 

For simplicity, consider only model (Q) with a single qubit as described in Sec. IVA. A bang-bang procedure applies 
a succession of resonant radiofrequeney pulscs [7r-pulses] which evolve one eigenstate of the qubit into the other, i.e., 
I t) — * I 4) an( l I ï) ~ * I T)> 011 a time scale Tp short compared to the bath correlation time. As before, throughout 
the following we assume the Schroedinger picture. If the rf field is strong enough, the interaction of the qubit with 
the environment can be neglected during the pulse, so that the bath coordinates remain unaffected. Therefore a total 
environment-entangled state, which reads [Eq.(p3|) with a convenient rearrangement] 

p{t-) - IT) /Wt (Tl + II) Pe,u (II + IT) Pe,n (II + II) ÍPe,n) 1 (Tl 
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immediately before the pulse, transforms into 

KU) = II) PeM UI + IT) Pe,ll (Tl + II) Pe,U (Tl + |T) (Pe.Tl) 1 (II 

immediately after the pulse. An elementary spin-flip cycle consists of two such pulses applied at a time interval Aí, 
the first of which reverses the qubit state, while the second restores the original configuration. Between any two pulses 
the qubit-environment system resumes the dynamics described by Hamiltonian (jï]). For the purpose of illustration, 
it is sufficient to examine the process in the limit of infmitely narrow pulses, Tp — ► |y, when each cycle can be 
approximated by the following piecewise evolution: i) propagation under Hamiltonian (Q) for a duration Aí, starting 
at time í„; ii) instantaneous 7r-pulse and interchange of qubit eigenstates states at time (t n + Aí); iii) evolution under 
Hamiltonian ([!]) from time (í„ + Aí) to time t n+ i — t n + 2Aí; iv) second 7r-pulse and interchange of qubit eigenstates 
at time í n +i. 

The analysis of this process in the framework of solution (^]) is quite straightforward provided we shift focus from 
the qubit eigenstates to the entangled environment distributions. Indeed, the transformation of a qubit-environment 
state upon application of a 7r-pulse can be interpreted also as a switch of the environment distributions associated 
with the qubit states, in the sense that /S e ,n * Pe n ) Pe,Tl ~^ Pe,Tli etc. The latter implies that the associated bath 
displacements before the rf pulse for Pe.TT' e ^ c -' become initial displaccments for p e ,ü, etc, after the pulse and vice 
versa, i.e., 

P u (nAt + 0) = (3 lq (nAt - 0) , 

and 



P lq (nAt + 0) = P u {nAt - 0) . 

For this reason, the interchange of qubit eigenstates following a 7r-pulse is equivalent to a rephasing of the entangled 
environment distributions. 

The time-dependence of the bath displacements for the entangled p e -s ensues now without difficulty. Starting at 
time t n — 2nAí [ío = 0] with displacements Pm) q (t n ), one has [see Eq.([26[)]: 

i) for t n < í < (í„ + Aí) = (2n + l)Aí: = [/%(T)g(*n) ± Xq/^q] exp[-iw g (í - t n )] =F X q /huj q ; 

ii) for (2n + l)Aí < í < t n+1 = 2(n + l)Aí: /3 T (|) ? (í) = [0i{r) q (t n + Aí) ± Xq/^q] exp[-íuj q (t -t n - Aí)] t Xq/^q- 
Furthcr, let the initial qubit-environment state be unentangled, with the environment in thermal equilibrium, such 
that /% 9 (0) = /3i q (0) — 0. Since the unperturbed evolution in the first half-cycle drives same mode displacements to 
vàlues of opposite sign, such that (i] q {t) = — /3j. g (í) = (x ? /^ 9 )[exp(— ioj q t) — 1] for < í < Aí, the corresponding 
evolved displacements will have opposite signs at any later time, that is, I3^ q {t) = — /3j. g (í) for all í. As a result, all 
phase factors [Eqs.(p5|a) and ([30|a)] in the bath correlation factor r)(t) for the reduced qubit state (Q) vanish, and the 
only relevant quantity remains the dissipative factor [Eq. ([30)3)] 

r(í) = 2^|/3 Tï (í)| 2 coth(?kV2fc b T) ■ (52) 

At the same time, the effect of a 7r-pulse is seen to amount to a change of sign of the bath displacements in the 
entangled environment distributions, such that 

P Ui)q {kAt + 0) = ~P ul)q (kAt - 0) . 

Accounting for this into the time-dependence of the displacements, yields immediately stroboscopic recurrence rela- 
tions of the form 

/3 u (kAt + 0) = u {{k - 2)Aí + 0)e^ UJ " At + ( e ^ l ^ At - l) 2 . (53) 

Tiu) q 

Solving the recurrence for k = 2n with initial condition /3| g (0) = recovers the stroboscopic result for T(í) obtained 
in ref . B , which we write in the form 



r sírob (2nAí) = 2^coth C^^j /3 T? , (2nAí) tan (^^J 



2 



where /3| g ,o(í) = (Xg/S-w 9 )[exp(— íuj q t) — 1] denotes the evolved displacement in the absence of pulses. Under the usual 
prescription for the density of states [see previous section] , the above expression becomes 
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- strob 



(2nAí) = Xh 2 Lj c J duj (j^j exp coth (gj^f) ^Tï.oC 2 ^*) tan (^^j 



(54) 



The decay of the corresponding bath correlation factor [usually referred to as qubit coherence] T] s trob 

(2nAí) = 

exp[— T st rob(2nAt)] was shown to be strongly suppressed at high enough pulse freqüències [w c Aí < 1], and to become 
completely quenched in the limit of continuous flipping. The effect is interpreted in ref. [p| as an approximate time 
reversal of the qubit evolution induced by each 7r-pulse, which becomes visible as an effective change of sign of the 
total Hamiltonian. The present point of view adds that the rephasing of bath entanglement responsible for the partial 
time reversal occurs precisely through a change of sign of the bath displacements and so amounts to a change of sign 
of the time derivative of the dissipative factor, since 



dT \ - / f\ií)„ 

—— = 2i > lü„ coth — — 7= 
dt ^ q \2k B T 



(55) 



In contrast, the exact time reversal of T(í) at time t rcquires that the bath displacements change according to 
y!/3| q) 9 (í+) = Xq[Pui)i^~)^* ■ Remarkably, this differs from the flip-induced change by a mere phase factor [see also 
M]. Thus if the bath correlation factor 77 = exp(— T) decays during the last moments of the first half-cycle, a change of 
sign of the displacements after the mid-cycle flip suffices to induce a subsequent revival of coherence. Formal evidence 
for this phenomenon is provided by the expression of the time-derivative dT / dt at the moment immcdiately following 
the mid-cycle pulse, which can be shown straightforwardly to amount to 

00 d— 1 

§((2» - l)M + 0) = -XKW Jdu (£) exp coth (^) cos 2 ((„ - i)o,Aí) tan (^) (56) 



for a Standard density of states. It is not difficult to observe that for d < 3 this expression is always negative, since the 
dominant contribution comes from the range < uAt < 7r, where the integrand is negative [note that for wAí > 7r 
the periodical factor cos 2 ((n — ^)wAí) tan {^^-) is multiplied by a monotonously decreasing function if d < 4]. Hence 
the correlation factor 77 always increases after the mid-cycle flip. 

Surprisingly, the detailed solution reveals that the second flip also induces a revival. Indeed, the derivative of T at 
the moment immediately preceding the second pulse of a cycle reads 



oc 



f (2nAí - 0) = Xh^ c J d. exp (-£) coth (^) sin^Aí) tan (^) (57) 



and, by the same argument as for expression (|5^), is noted to be always positive for d < 3, since the integrant is 
positive in the dominant range < cjAí < n. Because rj{t) is necessarily decreasing at this moment, the second 
pulse can only produce a revival. It also becomes apparent that rj displays [at least] a maximum between any two 
consecutive pulses. 

This effect is confirmed by a numerical integration of the exact dissipative factor ((5^), which becomes, for the 
selected density of states, 

T(í) = AftV / du (-) exp (--) coth (^=) \Mt)\ 2 . (58) 



uj c I V oj c J \2k B T 







The integration is straightforward, since the free evolution of the displacements P-\ q (t) is now known at all times 
through Eq.([2rj), while each rf-pulse contributes a simple change of sign in all displacements. Fig. [ijshows the result 
for representative model parameters, cycle periods and temperatures in an ohmic environment, d = 1 [the superohmic 
d = 3 case is qualitatively similar]. Should the time reversal of the bath correlation 7?(í) be exact, the màxima would 
oceur precisely at times t n = 2nAt, since the derivative of T is null in the initial state. Because the reversal is only 
partial, the màxima of rj(t) are seen to shift gradually, with each cycle, toward the midpoint of the interval between 
consecutive pulses. Moreover, as the temperature increases, the read-out vàlues at t n — 2nAt tend to become com- 
parable to the lowest vàlues in a cycle. One may infer that shifting the read-out times by ~ +Aí/2, so as to take 
advantage of these màxima, may yield a slight improvement of the overall outeome. 
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It turns out, in fact, that a minor rearrangement of the protocol results in a more significant gain. Consider the 
following version of the idealized bang-bang cycle: i) propagation under Hamiltonian (Q) for a duration Aí/2, starting 
at time t n \ ii) first 7r-pulse and interchange of qubit eigenstates states at time (í n + Aí/2); iii) evolution under 
Hamiltonian (Q) from time (t n + Aí/2) to time t n + 3Aí/2; iv) second 7r-pulse and interchange of qubit eigenstates 
at time t n + 3Aí/2; v) readout at t n +i = t n + 2At. This sequence is known in the literature on NMR decoupling 
methods as a symmetrized Carr-Purcell protocol [see ref. |Ï2| for the quantum computation context] , and was shown 
in the framework of the average Hamiltonian theory [AHT] to yicld an improvement in decoupling [suppression of 
decoherence] of order ~ 0(w c Aí) 4 over the Standard sequence. Since in the AHT this result is derived directly for 
the stroboscopic coherence [proportional to our bath correlation factor] at the end of each cycle, we find it instructive 
to review here a few additional details supplied by the exact coherent- pro duet solution to the associated relaxation 
problem. 

As for the Standard protocol, it can be verified that under this symmetrized sequence the qubit coherence also 
displays a maximum between any two consecutive rf-pulses. But since the read-out is scheduled now halfway between 
two pulses, one can expect that the corresponding output vàlues fall close to the coherence màxima following the 
second pulse in each cycle. In addition, we can anticipate that the symmetrized protocol also benefits from the halved 
lag [Aí/2] between the initial configuration to be preserved and the first applied rf-pulse, which induces the first 
revival. Similarly to an increase in the cycle frequeney, this reduced lag leads to faster clipping of the decoherence 
periods within each cycle, hence again to higher coherence màxima and higher read-out vàlues. 

For a quantitative assessment of these effects, let us resort once more to solution (|j). If the read-out times are 
t n = 2nAt, the bath displacements between consecutive 7r-pulses read as follows: 

i) for (í„ - Aí/2) < í < (í„ + Aí/2): /% )g (í) - - Aí/2 - 0) ± Xg/^g] exp[-iu q (í — í„ + Aí/2)] =p XqfhtjJq, 

ii) for (í„ + Aí/2) < í < í„ + 3AÍ/2: /3 ni)q (t) = + Aí/2 - 0) ± Xg/^g] exp[-ico q (í — t n — Aí/2)] =p Xg/^ 



As before, we assume an unentangled environment in thermal equilibrium, so that /3| ? (0) = /3j.ç(0) = and P\ q (t) 
-/3t q (t) for all í. Taking this into account leads to a stroboscopic recurrence of the form 

2ÍLj q At i ( -üjgAt/2 _ i\ 2 ( —iu q At 



1- 



P u (t n+1 ) = /3 Tï (í„)e- 2 ^ At + |j- (e"^ Aí / 2 - l) 



which solved for (3-\ q (0) = yiclds, under a Standard density of states, a stroboscopic dissipative factor 



2 f ( u> \ 1 f ui \ ( hui 
Y sym (2nAt) — Xh lü c I duo — ) exp ) coth 

o 



2k B T 



(3 u , Q (2nAt) 



1 - cos(w 9 Aí/2) 



cos(w 9 Aí/2) 



2 



(59) 



Here /3f 9i o(2nAí) is again the unperturbed evolved displacement. Note that expression ( p9[ ) differs from the 
Standard expression only through the substitution of the factor | tan(w g Aí/2)| 2 by a factor of |[1 — 

cos(oj g Aí/2)]/ cos(uj q Aí/2) | 2 . Remarkably, the two factors have a qualitatively similar functional dependence on 
(ujqAt), although their periodicity differs from 2tt to 47r, respectively. Quantitatively, on the dominant range 
< uiqAt < 7r, one has | tan(w 9 Aí/2)| 2 > |[1 - cos(w g Aí/2)]/ cos(w 9 Aí/2)| 2 and this proves sufHcient to render 
r s íro6(2n.Aí) > T sym {2nAt), hence rjstrobi^nAt) < rj sym (2nAt). The magnitude of this effect has been evaluated 
by numerical integration of the stroboscopic expression (^9|) and the results are displayed in Fig. ^, alongside the 
corresponding Standard output. In addition, a numerical integration of expression ( p8|) for the corresponding exact 
displacements at intermediate times shows clearly [Fig. |^] the emergence of the improved suppression of decoherence 
from the actual revival of coherence induced by each rf-pulse. 

As expected, the qubit coherence is always better preserved under the symmetrized protocol. In addition, the 
symmetrized outeome is considerably more stable under temperature changes, with coherence variations < 10~ 4 
compared to ~ 10~ 2 under the Standard protocol. The absolute magnitude of the improvement generally does not 
exceed 5% [over time] at identical cycle periods and temperatures in an ohmic environment, but may raise above 
~ 10% in the superohmic case [not shown]. However, the true advantage of the symmetrized sequence can be better 
appreciated in the high precision regime. Indeed, as can be observed from Fig. ^, in order to maintain the errors in the 
bath correlation r\ under, e.g., ~ 0.1% over an extended period of time [say, at least an order of magnitude longer than 
the typical decoherence time ~ w" 1 ], a Standard sequence must be applied at a frequeney substantially higher than 
u) c . To give a semiquantitative reference, under the choice of model parameters employed in Fig. ^ the rf-frequeney 
must be at least 4w c at T = 1K and at least 6ui c at T = Í00K. The symmetrized protocol is seen to deliver a 
similar precision with a sizable reduction in frequeney, in an essentially temperature-independent manner. E.g., under 
the same model parameters the necessary frequeney amounts to only ~ 2ui c at both T = 1K and T — 10QK . This 
reduction in the operating frequeney increases as the precision bounds on the correlation rj^ increase. As a result, 
a symmetrized protocol applied at T = 100-R" with a frequeney of 2.5 ui c can perform better [errors < 10~ 4 ] than a 
Standard protocol operating at T = IK with a frequeney of 10w c [errors < 10~ 3 ]. 
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VII. CONCLUSION 



In principle, the dynamics of a quantum register on a direct sum of DFSs may remain decoherence free even 
when the reduced register state is not restricted to a single DFS and involves considerable entanglement with the 
environmcnt. The present paper provides explicit conditions for the occurrence of this effect in a large class of thermal 



coherent-product distributions [Eq.(24)], supported by quantum registers with DFSs under interactions linear in the 
coordinates of a harmònic environment [Eq. @]. Such coherent-product states represent distributions on direct 
sums of register DFS, nontrivially entangled at all times with a statistical superposition of environmental Gaussian 
states. In particular, any register distribution on a direct sum of DFS develops precisely into a coherent-product 
distribution when brought in contact with an initially uncorrelated thermal environment. On the other hand, in the 
zero-temperature, pure state limit these states become self-consistent [but not soliton bearing] Davydov coherent- 
product states. 

Our study was prompted by the straightforward observation that the total state described by Eq.fpij) yields a 
unitary register evolution in the special case when the associated bath displacements coincide with the stationary dis- 
placements characteristic of each register DFS. Here we showed that the set of coherent-product states that propagate 
the reduced register state in an unperturbed, decoherence-free manner is considerably wider. These states require an 
appropriate density of degenerate environmental modes such that i) the displaced harmònic environments experienced 
within contributing DFSs have identical energy reference points [Eqs.(|36|a) or (JÏ5|a)]; ii) the environment Gaussian 
distributions can be adequately synchronized through their Gaussian displacements [E qs. (|37|)] ; and iii) the amplitude 
of the environment mediated correlations between diffcrcnt DFS is finite [expression (g^) must be finite]. The latter 
requirement implies that the displaced environment vacua seen by contributing DFS must be unitarily equivalent. 
Somewhat unexpected, an ohmic environment cannot support such states in single qubits or registers with weak 
collective decoherence, but is marginally effective for linear registers with individual decoherence. A superohmic en- 
vironment, on the other hand, is considerably more efficicnt in sustaining entangled, decoherence-free distributions in 
a variety of registers, including registers with weak collective decoherence. However, since the amplitude of correla- 
tions between distinct DFS compatible with such a propagation decreases significantly with increasing temperature, 
nontrivial entangled decoherence-free states likely entail a low temperature environment. The simplest examples of 
joint register-environment states supporting multi-DFS decoherence free register propagation are Davydov ansatz 
pure states with properly phased coherent environment components. 

The practical worth of these distributions depends to a large extent on one's ability to generate and control a proper 
phasing of the environmental entanglement. It turns out that qubit flipping through radiofrequeney 7r-pulses, as in 
bang-bang suppression of decoherence, provides an elementary procedure for such a manipulation, which amounts to 
a change of sign of the associated Gaussian displacements. As a corollary, examination of the bang-bang cycle in 
terms of the coherent-product states (|24|) reveals that each rf-pulse in the protocol causes a revival of coherence. It 
also follows that the well-known time-symmetric version of the cycle takes advantage of these revivals in a natural 
way, by approximately matching the read-out times to the times of maximum qubit coherence between consecutive 
pulses. As a result, high precision control becomes possible with cycle freqüències reduced by a factor of about 2-^4 
and with a temperature sensitivity diminished by two orders of magnitude. 

Returning to DF coherent-product distributions, let us assume the availability of a realistic preparation procedure. 
The coherent concatenation of distinct register DFSs effected by such states immediately suggests the possibility of 
expanding the capacity for noiseless storage beyond the limit set by parallel [incoherent] storage on the same set of 
DFSs. Consider, as a simple example, a 3-qubit register with weak collective decoherence as described in Sec.VB. The 
register Hilbert space decomposes into two 1-dimensional DFS, corresponding to each of the states | |||) and | Tíí), 
and two 3-dimensional DFSs spanned by the basis sets {| J,jt)> I iíi) I TU)} an d {| Iíí), | TIT) | ííi)}, rcspectively. 
According to the usual theory of DF propagation, each of the 3-dimensional DFSs supports noise-protected storage 
of one encoded qubit or qutrit. In contrast, a coherent-product state, e.g., a Davydov state of the type described in 
Sec.IIIB, reading 



KitI Uí> + ami iíi) + «tul TU)] ® l/M*)) + KttI iíí) + «titI íií) + «ml TU)] ® |A(*)) 



= ai|tfi}<8>|/?i(i))+a2|*2)® |/? 2 (i)> , 

allows the noiseless storage of two entangled encoded qubits or qutrits [that is, of a genuine two-qubit or two-qutrit 
memory] when the coherent environment states |/?i i2 (í)) have a finite, time-independent overlap, (/3i(í)|/?2(í)) = 
const. ^ 0. Notably, the latter condition makes it possible to extract the relative phases of states belonging to 
diffcrcnt DFSs [li'l, 2) above] solely by means of register measurements, provided the overlap {(3i\(Í2) can be determined 
indepcndently. For instance, one may use a test state of the form 
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i*> = -^[iut>®iA(í)> + im>® !/&(*)>] , 

and projective measurements on the register states |</>) = (1/ a/2) [| I1T) + I ITT)] an( l \4>') = (l/\/2)[| TIT) + A TTT)] 
to extract the desired overlap from the corresponding probabilities p<p = (cf)\TrB[\&)(&\]\<j>) = 1/2 + Re(j3i\(32) and 
P0/ = ((/)'|rrB[|$)($|]|0') = 1/2 + /m(/3i|/3 2 ), respectively. Similar register measurements on any encoding state 
l'S) will then allow a full read-out of the information stored in the amplitudes ajki- As usual for a mixed register 
state, a complete set of such measurements can be obtained from, e.g., projective measurements on each of the states 
in the computational basis, plus measurements on all entangled superpositions of pairs of basis states of the form 
(l/\/2)[|a) + and (l/\/2)[|a) + i\(3)]. In general, as shown in Sec.VB, noiseless coherent-product superpositions 
of DFS under weak collective decoherence enable a doubling in size of the code space, corresponding to the addition 
of an extra encoded qubit. Note however that the advantage of coherent-product encoding becomes fully apparent 
on mixed distributions, where the set of independent storage liberties [matrix elements] becomes considerably larger 
than in the pure state case. 

We conclude by pointing out that the environment-entangled unperturbed evolution discussed in this paper should be 
understood as one particular case of environment-entangled wiiíarj/ dynamics. For instance, conditions (p6p) or system 
( |37| ) sufEce to guarantee a unitary, quasi-unperturbed propagation of pn (t) . This is because the stationary phase term 
Fq in Eq. ( |33| ) is separable into contributions from individual, orthogonal DFS, and contributes merely stationary 
energy shifts to the unperturbed Hamiltonian. Furthcrmore, it is also possible to obtain sufficicnt conditions for the 
unitary propagation of pFt{t) by requiring only that the dissipative factors of all ?7j!/,jv be stationary [alT j v ,j< v < / dt = 0] 
and that all phase factors be separable, which amounts in fact to d$ v / 'dt = 0. The unperturbed Hamiltonian is 
modified then by time- dependent energy shifts and the bath displacements must satisfy the constraints 

E (tfU - J (ti, - ti>«) = * 
i 

and 

E - = o . 

The requirement that the stationary dissipative factor Tj j, be finite remains, of course, unchanged. 
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Figure captions 



FIG. 1. Time-dependence of qubit decoherence in an Ohmic environment under Standard [hollow symbols] and sym- 
metrized [sòlid symbols] sequences of rf-pulses, for two typical temperatures corresponding to ksT /Tujj c — 0.01 [squares] and 
UbT/UíOc = 1.0 [circles]. For a cut-off frequency lo c ps 10 13 s _1 , the respective temperatures are T = 1K and T = 100iV Time 
is given in units of the read-out [cycle] period r cyc i e — 2Aí, and the model scaling constant was set to Xh 2 uj c = 0.25. 



FIG. 2. Errors in qubit decoherence as a function of cycle frequency under Standard [Eq. (p4|), hollow symbols] and sym- 
metrized [Eq. (^), sòlid symbols] sequences applied for a time t, at temperatures T = 1K [squares] and T = W0K [circles]. 
Other parameters as in Fig. Q Each point represents the read-out after a number of t/T cyc i e cycles. 
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